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We construct new models of black hole-neutron star binaries in quasiequilibrium circular orbits 
by solving Einstein's constraint equations in the conformal thin-sandwich decomposition together 
with the relativistic equations of hydrostationary equilibrium. We adopt maximal slicing, assume 
spatial conformal flatness, and impose equilibrium boundary conditions on an excision surface (i.e., 
the apparent horizon) to model the black hole. In our previous treatment we adopted a "leading- 
order" approximation for a parameter related to the black-hole spin in these boundary conditions 
to construct approximately nonspinning black holes. Here we improve on the models by computing 
the black hole's quasilocal spin angular momentum and setting it to zero. As before, we adopt 
a polytropic equation of state with adiabatic index F = 2 and assume the neutron star to be 
irrotational. In addition to recomputing several sequences for comparison with our earlier results, 
we study a wider range of neutron star masses and binary mass ratios. To locate the innermost 
stable circular orbit we search for turning points along both the binding energy and total angular 
momentum curves for these sequences. Unlike for our previous approximate boundary condition, 
these two minima now coincide. We also identify the formation of cusps on the neutron star surface, 
indicating the onset of tidal disruption. Comparing these two critical binary separations for different 
mass ratios and neutron star compactions we distinguish those regions that will lead to a tidal 
disruption of the neutron star from those that will result in the plunge into the black hole of a 
neutron star more or less intact, albeit distorted by tidal forces. 

PACS numbers: 04.30.Db, 04.25. Dm, 04.40.Dg 



I. INTRODUCTION 

Coalescing black hole-neutron star (BHNS) binaries, 
as well as other compact binaries composed of neutron 
stars and/or black holes, are among the most promis- 
ing sources of gravitational waves for both ground-based 
P) S) H) 0] and space-based laser interferometers [1, [1] . 
BHNS binary mergers are also candidate central en- 
gines of short-hard gamma-ray bursts (SGRBs) (see, 
e.g. 01 and references cited therein). The remnants 
of both BHNS binary merg ers Q and binary neu- 
tron star mergers [Tol. I 111. [l3 . [Tsj are feasible progenitors 
for SGRBs because both may result in black holes sur- 
rounded by hot, massive accretion disks with very little, 
if any, baryon contamination along the polar symmetry 
axis. 

Motivated by these factors, considerable effort has 
gone into the study of BHNS binaries. Most approaches 
to date assume Newtonian gravity in either some or all 

tects of the calculation (see, e.g. [H, [l^ [H, |T3, 
[20I [21 ^2 ^ 23I f or quasiequilibrium calculations and 
[2j,l2in2l~|23rii, Ei, H in for dynamical simula- 



tions). More recently, several groups have also studied 
BHNS binaries in a fully relativistic framework, both for 
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quasiequilibrium models |33. l33l. |3J. [33, [30, [33| and dy- 
namical simulations [sl. [ol l38l l39l. |40||. 

Our group has pursued a systematic approach to de- 
veloping increasingly realistic models of BHNS binaries 
in quasiequilibrium circular orbits. Our first studies 
d, [H, [H, [11] assumed extreme mass ratios, i.e., black 
hole masses that are much greater than the neutron star 
mass. While this is a very natural first step from a 
computational point of view, binaries with comparable 
masses are much more interesting from the perspective 
of ground-based gravitational wave observations and for 
the launching of SGRBs. More recently we have therefore 
relaxed this assumption and have extended our results to 
the case of comparable- mass BHNS binaries [H, [13] ■ 

Specifically, in [U (hereafter Paper I) we constructed 
quasiequilibrium models by solving Einstein's constraint 
equations in the conformal thin-sandwich formalism, as- 
suming conformal flatness and maximal slicing, together 
with the relativistic equations of hydrostationary equi- 
librium. We accounted for the black hole by excising a 
coordinate sphere and imposing the equilibrium black- 
hole boundary conditions of Cook and Pfeiffer |41j] . This 
original version implemented a "leading-order" approx- 
imation to nonspinning black holes, which equates an 
otherwise undetermined spin parameter f2j. that appears 
in the boundary condition for the shift vector with the 
orbital angular velocity seen by an incrtial observer at in- 
finity, f2. As for the original irrotational binary black hole 
models of [4l|, this condition does not lead to simulta- 
neous turning points of the binding energy and the total 
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angular momentum in constant-mass sequences in Pa- 
per I. Such simultaneous turning points are expected for 
those sequences if they are truly in quasiequilibrium [13] . 
An improvement over this condition, namely to iterate 
over Qr until the quasilocal spin angular momentum of 
the black hole vanishes, was suggested and implemented 
for binary black holes by [1^. 

In this paper we reconstruct quasiequilibrium models 
of BHNS binaries using the same techniques as in Paper 
I, but with the improved black hole spin angular veloc- 
ity condition as suggested by [11]. We then compute 
sequences of BHNS binaries in quasicircular orbits for 
a wider range of neutron star masses and binary mass 
ratios than in Paper I, focusing our attention on irro- 
tational neutron stars orbiting nonspinning black holes. 
Here we focus only on the irrotational state for the neu- 
tron star because it is astrophysically considered to be 
more realistic in a BHNS binary [ll, [H, [11]. On the 
other hand, we will compute the case of spinning black 
holes in future work. As was the case for the irrotational 
black hole binaries constructed in [i^, we find that this 
improved condition for the spin parameter of the black 
hole fir does lead to simultaneous turning points in the 
binding energy and the total angular momentum along 
constant-mass sequences. 

The paper is organized as follows. We briefly review 
the basic equations in Section II. We present numerical 
results in Section III, and outline some qualitative consid- 
erations concerning the fate of BHNS binaries in Section 
IV. In Section V we summarize our findings. Throughout 
this paper we adopt geometrized units with G = c ~ 1, 
where G denotes the gravitational constant and c the 
speed of light. Latin and Greek indices denote purely 
spatial and spacetime components, respectively. 

II. FORMULATION 

In this Section we review the equations we solve to con- 
struct a quasiequilibrium BHNS binary. The equations 
are very similar to those in Paper I, but, as we explain 
below, we have improved both our algorithm for the solu- 
tion of the gravitational field equations and the method 
to compute the spin parameter fir that appears in the 
black hole boundary condition for the shift vector. For 
a more detailed discussion about the formalism we use 
here, we refer to Paper I, the review articles [4l,[43|, and 
Sec. II of [1^ for the hydrostatics. 

A. Gravitational field equations 

The line element in 3 + 1 form is written as 

= -a^dt"^ +-f^j{dx' + (3'dt){dx^ + /3^dt), (1) 

where g^^, is the spacetime metric, a the lapse function, 
the shift vector, and jij the spatial metric induced 



on a spatial slice S. The spatial metric jij is further 
decomposed according to 7^ = ip'^^ij, where ip denotes 
the conformal factor and jij the conformal background 
spatial metric, defined such that det 7 = 1 in Cartesian 
coordinates. We also decompose the extrinsic curvature 
AT'^ into a trace K and a traceless part A'^^ according to 

K'^ = ip-^^A'^ + \-i'^K. (2) 

The Hamiltonian constraint then becomes 

= -2^VV + \i^R + i^V^'A-' - i^-'^i^^. (3) 

Here = 7*-' ViVj is the covariant Laplace operator, 
Vi the covariant derivative, Rij the Ricci tensor, and 
R — Rij the scalar curvature, all associated with the 
conformal background metric 7^ . 

We employ the conformal thin-sandwich decomposi- 
tion of the Einstein equations [1^. In this decomposi- 
tion, we use the evolution equation for the spatial metric 
to express the traceless part of the extrinsic curvature in 
terms of the time derivative of the background metric, 
iiij = dt"/ij, and the gradients of the shift vector. Under 
the assumption of equilibrium, i.e., Uij = in a corotat- 
ing coordinate system, the traceless part of the extrinsic 
curvature reduces to 

iy = f. f^^pj + yjp^ _ If^S/.p'^) . (4) 
2a V 3 / 

Inserting Eq. Q into the momentum constraint we ob- 
tain 

V^P' + lv'{\7jf3^) + R)(3^ 
3 

= 167raV'*/ -t- 2A'^Vj (a^^^) + ^a7*^ AT. (5) 

For the construction of quasiequilibrium data it is also 
reasonable to assume 9tAr = in a corotating coordi- 
nate system. The trace of the evolution equation for the 
extrinsic curvature then yields 

+ai}~^AijA'^ - 2f^V,aVj In (6) 

The matter terms on the right-hand side of Eqs. jSj, 
([5|), and ([6]) are derived from the projections of the stress- 
energy tensor T^;^ into the spatial slice E. Denoting the 
future-oriented unit normal to S as n^, the relevant pro- 
jections of T^^ are 



p = n^n,T^'\ (7) 

f = ~7>.^^^ (8) 

Srj - l^^ljuT^", (9) 

S = f''S,,. (10) 
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Here we write the stress-energy tensor as 



(po + Pi + P)Uf_,u,y + Pg^,^, 



(11) 



assuming an ideal fluid. The quantity is the fluid 4- 
velocity, po the baryon rest-mass density, pi the internal 
energy density, and P the pressure. 

Equations (O, (O and ^ provide equations for the 
lapse function a, the shift vector /?*, and the conformal 
factor ip, while A^^ can be found from Eq. The con- 
formally related spatial metric 7y and the trace of the ex- 
trinsic curvature K remain freely specifiable, and have to 
be chosen before we can solve the above equations (note 
that we have already set the time derivatives of these 
quantities, which are also freely specifiable, to zero). As 
in Paper I we assume a flat background jij = rjij , where 
Tjij denotes a flat spatial metric, and maximal slicing, 
K — 0. In Cartesian coordinates, Eqs. ([3]), ([5]) and © 
can then be written as 



(12) 



+ -d'idjP^) = IGn^ij^f + 2i*^'9j($V"^)(13) 

A$ = 27r$^'*(p -I- 2S) + -<^^/r^A,jA'3, (14) 

8 

where A and di denote the flat Laplace operator and the 
partial derivative, and $ = aip. Eq. ([4]) becomes 



(15) 



For numerical purposes we further decompose the vari- 
ables and their equations into parts associated with the 
black hole and the neutron star. For details of this de- 
composition we refer to Appendix A of Paper I fs?'] . 

In Paper I we solved directly for the lapse function a, 
the conformal factor and the shift vector Instead 
of solving Eq. ^ for the lapse, we now solve Eq. (|14p for 
the combination $ = aip. This choice is quite common 
(see e.g. [i^ 113), and has the advantage of eliminating 
the source term —2rfWiadjh\il} on the right-hand side 
of Eq. ^ above, or equivalently Eq. (9) in Paper I. 
This term falls off like 1/r'* , and hence more slowly than 
AijA^^ . Eliminating this term therefore enhances the ac- 
curacy of the numerical solution. We will quantify the 
improvement in Section [III Bl 



B. Boundary conditions 

In order to solve the gravitational field equations p2| , 
(|13p . and (jl4p . we have to impose appropriate boundary 
conditions on two different boundaries: outer boundaries 
at spatial infinity and inner boundaries on the black hole 
horizons. 

The boundary conditions at spatial infinity follow from 
the assumption of asymptotic fiatness. With the help of a 
radial coordinate transformation u = 1/r in the external 



computational domain, our computational grid extends 
to spatial infinity [isl . [soj , and we can impose the exact 
boundary conditions 



= 1, 

, = {^x R)\ 
= 1. 



(16) 
(17) 
(18) 



Here is the orbital angular velocity of the binary system 
measured at infinity, and R = (X, Y, Z) is a Cartesian 
coordinate centered on the center of mass of the binary 
system. We express the shift vector in a corotating 
coordinate system that we adopt throughout our calcu- 
lation. In an inertial coordinate system, the shift vector 
would tend to zero at spatial infinity, while in the coro- 
tating coordinate system of the numerical code the shift 
vector diverges at spatial infinity. For computational pur- 
poses, it is therefore convenient to write the shift vector 
as a sum of the rotational shift term j3l^^ = (il x R)"^ and 
a residual part (which tends to zero at spatial infinity), 
and solve the equations only for the latter. 

The inner boundary conditions arise from the excision 
of the black hole interior. The assumption that the black 
hole is in equilibrium leads to a set of boundary condi- 
tions for the conformal factor and shift vector [4l| (see 
also [H, [m as well as the related isolated horizon for- 
malism, e.g. [H, H^). The boundary condition for the 
conformal factor is 



5 4 '' 



(19) 



where = "0 is the outward pointing unit vector 
normal to the excision surface and /ly is the induced 
metric on the excision surface, hij = ip'^hij = 7ij — SiSj. 
The quantity J is computed from the projection of the 
extrinsic curvature Kij as J = h'^Kij. The boundary 
condition on the normal component of the shift vector is 



P±\s = a\s- 



(20) 



The tangential components must form a conformal 
Killing vector of the conformal metric hij on the exci- 
sion surface (see [4l|). This can be achieved by choosing 
them to be Killing vectors of a 2-sphere, 



"ok 



(21) 



Here ftl is a freely specifiable vector, related to the black- 
hole spin, that we take to be aligned with the Z-axis, and 
x'^ is a Cartesian coordinate centered on the 2-sphere. 

We assume the excision surface to be a coordinate 
sphere. In our previous treatment we implemented 
a "leading-order" approximation to nonspinning black 
holes and set il^. = fi, where is the orbital angular 
speed (compare JUi] ) . Following ^4^ we now iterate over 
rir until the black hole's quasilocal spin angular momen- 
tum 



(22) 
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vanishes. Here is an approximate Kilhng vector of 
hij that we find by solving the Killing transport equa- 
tions as described in [i^l (see also [H^ for more detailed 
descriptions and [s&l for alternative methods for finding 
approximate Killing vectors of closed 2-surfaccs). 

According to ^3] , the boundary condition on the lapse 
function can be chosen freely. In this paper, we choose a 
Neumann boundary condition 







(23) 



on the excised surface. 

We refer to Sections H.E and H.F of Paper I for a 
discussion of how the orbital angular velocity, the center 
of rotation, and several global quantities including total 
angular momentum and mass are computed. 



C. Numerical Method 

As in Paper I [37j . we construct our numerical code 
based on the Lorene spectral-methods library routines 
developed by the Meudon relativity group [56|. In our 
code, the computational space is broken into multiple 
domains. Each domain around the black hole is covered 
by Nr X Ne X ^ Al X 33 X 32 or 49 x 37 x 36 col- 
location points, while those around the neutron star are 
covered by 25 x 17 x 16 collocation points. Here iV^, Ng, 
and A''^ denote the number of collocation points for the 
radial, polar, and azimuthal directions, respectively. We 
use a larger number of collocation points for the black 
hole domains than for the neutron star domains because 
the source terms of the black hole equations are sensitive 
to the resolution of the neutron star. Since the black 
hole domains are centered on the black hole, we need 
a higher angular resolution to adequately resolve these 
source terms. The neutron star equations, on the other 
hand, have large source terms only near the neutron star. 
Since the neutron star domains are centered on the neu- 
tron star, a more modest angular resolution is sufficient 
to resolve these terms. 

We refer to Appendix A of Paper I for a detailed dis- 
cussion of the decomposition of the equations and their 
source terms. 



III. NUMERICAL RESULTS 

Throughout this paper, we model the neutron-star 
equation of state by the polytropic relation P = kPq, 
where P is the pressure, po the baryon rest-mass den- 
sity, r the adiabatic index, and k a constant. We choose 
F = 2 for the adiabatic index, and compute several dif- 
ferent constant-mass inspiral sequences. Specifically the 
rest mass of the neutron star and irreducible mass of 
the black hole are kept constant along each sequence. 
We focus on baryon rest masses for neutron stars in the 
range of 0.12 < M^^ < 0.17 (see Fig. [J). Here and in 



the following we normalize all dimensional quantities in 
terms of the polytropic length scale i?poiy = k^^'--^^^^\ 
e.g. Mg^ = A/g^/iJpoiy. In terms of compaction, our 
models are in the range 0.1088 < C < 0.1780, where 



C = 



^"ADM.O 
i?0 



(24) 



-''-'adm.o 



is the compaction of a spherical neutron star 
the neutron star's ADM mass in isolation, and Rq its 
areal radius. Our most compact polytropic model with 
Afg^ — 0.17 is very close to the maximum baryon rest 



mass for spherical F 
0.180. 
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FIG. 1: Mass-radius relation for a spherical, polytropic star 
with adiabatic index F = 2. The horizontal axis denotes the 
isotropic radius and the vertical axis the baryon rest mass 
(solid curve) and gravitational mass (dashed curve), where the 
radius and masses are in polytropic units. The gravitational 
mass is equivalent to the ADM mass for an isolated, spherical 
neutron star. Filled circles denote the models we choose in 
this paper. The compaction parameter C for each model is 
also shown. 

We consider mass ratios in the range 1 < 9 < 10, where 
we define the mass ratio as 

^^J^^ (25) 

i.e., the ratio of the irreducible mass of the black hole 
(M^^) to the ADM mass of a spherical, isolated neu- 
tron star (M^^j^g). Here the irreducible mass is 
identical to the ADM mass for an isolated nonspinning 
black hole. Note again that we fix the irreducible mass 
of the black hole and the baryon rest mass of the neu- 
tron star for the construction of constant-mass sequences. 
For the definition of the mass ratio, however, we use the 
ADM mass of a spherical isolated neutron star Af^DM O' 
because this turns out to be more convenient for com- 
parisons with third-order post-Newtonian (3PN) results 

We tabulate our numerical results in Appendix [X\ 
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A. Configurations 
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FIG. 2: Contours of the conformal factor in the equatorial 
plane for our innermost configuration with binary mass ratio 
q — 3 and neutron-star baryon rest mass M^^ = 0.15. The 
cross "x" indicates the position of the rotation axis. 

We show contours of the conformal factor tjj for a 
BHNS binary with mass ratio q = 3 and neutron-star 
baryon rest mass Afg^ — 0.15 in Fig. [21 This figure 
represents the configuration at the smallest orbital sep- 
aration for which our code converged. The thick sold 
circle on the left-hand side denotes the position of the 
excised surface (the apparent horizon) , while that on the 
right-hand side denotes the position of the neutron-star 
surface. The value of on the excised surface is not 
constant because a Neumann boundary condition (jl9p is 
applied. 



The solid curve shows the new results in this paper, the 
dashed curve the results of the old formulation in Paper 
I but including the new method to compute Qr, and the 
dotted-dashed curve the old results in Paper I. For closer 
configurations (larger f2A/o), the difference is dominated 
by the accuracy of the spin parameter fir, while for larger 
separations (smaller flMo), the difference comes from the 
change of variables. 
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FIG. 3: Binding energy E^^/Mq (upper panel) and total an- 
gular momentum J/Mq (lower panel) as a function of Q,Mq 
for binaries of mass ratio q = 1 and neutron-star mass 
— 0.15. The solid line with filled circles shows new 
results, and the dashed line the old results in Paper I. The 
dotted-dashed line denotes the results of the 3PN approxima- 
tion [s^l • The numerical sequences end due to cusp formation 
- and hence the onset of tidal disruption - before the bi- 
nary reaches the innermost stable circular orbit at the turn- 
ing point of the binding energy. 3PN sequences all exhibit 
turning points and cannot reveal cusps. 



B. Binding energy and total angular momentum 

In Figs. [3] ~ [5] we show the binding energy Eh/Mo 
and the total angular momentum J/Mq as a function of 
the orbital angular velocity ^IMq for neutron stars with 
baryon rest mass Mg^ = 0.15 and for mass ratios g = 1, 
3, and 5. Here we define the binding energy as the dif- 
ference in the total ADM mass of the binary from that 
at infinite orbital separation, Ey^ = Madm — Mq, and the 
total ADM mass of the binary at infinite orbital separa- 
tion as Mq = + M^^y^ Q. In these figures we also 
include, for comparison, results from Paper I and 3PN 
approximations [57|. We find that our new results show 
much better agreement with 3PN results especially for 
larger separations (smaller flMg). This improvement is 
due to the change of variables we discussed in Section 
III Al In Fig. [6] we show the relative difference of the total 
angular momentum from that of 3PN approximation as a 
function of the orbital angular velocity flMo for neutron- 
star baryon rest mass Mg^ = 0.15 and mass ratio q = 5. 
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FIG. 4: Same as Fig. O but for sequences of mass ratio q — 3. 
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FIG. 5: Same as Fig. |3]but for sequences of mass ratio q = 5. 
The binary encounters an ISCO before the neutron star is 
tidally disrupted. 



FIG. 7: The binding energy as a function of total angular 
momentum for binaries of mass ratio q = 5, and different 
neutron star compactions. 
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FIG. 6: The relative difference between the 3PN and numeri- 
cal total angular momentum as a function of QMq for binaries 
of mass ratio q — 5 and neutron-star mass Mg^ = 0.15. The 
solid curve shows new results in this paper, the dashed curve 
the results of the old formulation in Paper I but using the new 
method to compute fir, and the dotted-dashed curve the old 
results in Paper I. This figure demonstrates that the new spin 
condition improves the results predominantly at small binary 
separation, while the new formulation of the equations has a 
larger efi'ect at large binary separations. 



Before discussing these results in greater detail it is 
useful to anticipate some qualitative scaling. From a very 
crude Newtonian argument, we can estimate the binary 
separation c?tid at which the neutron star will be tidally 
disrupted by comparing the tidal force exerted by the 
black hole on a test mass m on the neutron star's surface 
with the gravitational force exerted by the neutron star 



Jtid 



-^^NS \ ^/'^ RnS 



Ah 



BH 



NS 



(26) 



If dtid is greater (by a sufficient amount) than the inner- 
most stable circular orbit (ISCO) separation disco we 
may expect the neutron star to be tidally disrupted be- 
fore being swallowed more or less intact, albeit distorted 
by tidal forces, by the black hole. The qualitative rela- 
tion ([26]) suggests that the tidal separation decreases with 
increasing mass ratio q and neutron star compaction. 

Our sequences terminate shortly before the onset of 
tidal disruption. We therefore expect to encounter min- 
ima in the binding energy and angular momentum, which 
identify the ISCO, only for binaries with sufficiently large 
mass ratio q and neutron star compaction. Comparing 
Figs. [3] - [5] we indeed find turning points only for the 
largest mass ratio g = 5. 

We also note that in this case the turning points in the 
binding energy and the angular momentum occur simul- 
taneously to within our numerical accuracy. This was not 
the case with for our earlier results, which adopted the 
"leading-order" nonspinning condition fi,. = fJ. When 
imposing the more accurate condition S — Q for irro- 
tational binaries, we now do find simultaneous minima 
(compare also the analogous results of [4l|, for binary 
black holes). 

To highlight this finding we graph in Fig. [7] the binding 
energy versus angular momentum for sequences of mass 
ratio g = 5, but different neutron star compactions. A 
simultaneous turning point in the binding energy and an- 
gular momentum leads to a cusp in these curves. As sug- 
gested by Eq. (j26p , we do not find turning points for small 
compactions (since the sequences end at tidal disruption 
before encountering an ISCO), but for larger compactions 
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these curves indeed form a cusp. While our results agree 
with 3PN results very well, we do note a small deviation 
that increases with the neutron star compaction, as one 
might possibly expect. We also clearly find that our new 
data agrees with 3PN results much better than the old 
ones in Paper 1 for smaller binary separation (see curves 
of = 0.15 in Fig. [T]). We note that 3PN sequences 
cannot identify tidal disruption and therefore always ex- 
hibit turning points. 



C. Quasilocal spin angular momentum of the black 

hole 

Probably the most important numerical improvement 
that we present here is the incorporation of a method to 
compute the spin angular velocity of the black hole into 
our numerical code. We first obtain the Killing vector 
on the excised surface by solving the Killing transport 
equations, and then compute the quasilocal spin angular 
momentum psl [53 | . Requiring the angular momentum 
to be zero {S = 0), we iterate over the black hole spin 
parameter fir- In Paper I, by contrast, we simply set 
this parameter equal to the orbital angular velocity 17, 
resulting in a "leading-order" approximation. In Figs. 
IS] and [21 we compare the quasilocal spin angular mo- 
mentum and black hole spin parameter Qr, calculated 
under the two conditions fir = and 5* = 0. For both 
computations, we use our set of improved variables (as 
opposed to reusing our results from Paper I). In these 
figures, solid and dotted-dashed lines represent configu- 
rations with zero quasilocal spin angular momentum for 
the black hole, while dashed and dotted curves represent 
configurations where we set the spin angular velocity of 
the black hole equal to the orbital angular velocity. 

We find from Fig. [5] that the quasilocal spin angular 
momentum of the black hole was negative and increased 
in magnitude as the orbital separation decreased when 
we used the condition fi^ = li in our boundary condi- 
tions. This negative spin angular momentum results in 
a decrease in the total angular momentum compared to 
the S = case. This explains why we did not find a min- 
imum in the total angular momentum along a sequence 
in Paper I. Here, by maintaining S = within numerical 
errors, we find the minimum of the total angular momen- 
tum is coincident with that of the binding energy. 

Similarly, the fact that the quasilocal spin angular mo- 
mentum becomes negative for the condition Qr — ^ im- 
plies that the spin angular velocity fir should be smaller 
than the orbital angular velocity when we require 
5 = 0, as confirmed by Fig. O 



IV. QUALITATIVE CONSIDERATIONS 

Among the most important questions in the context of 
BHNS binaries is whether the coalescence leads to a tidal 
disruption of the neutron star, or whether the neutron 



Mg =0.15 



0- =7^^ 



-0.01 - 



-0.02 - 



-0.03 



-0.04 



0.02 



0.04 0.06 



0.08 



0.1 



FIG. 8: Quasilocal spin angular momentum of the black hole 
as a function of the orbital angular velocity (along a quasiequi- 
librium constant-mass sequence). Solid and dotted-dashed 
lines correspond to the condition S — 0, while dashed and 
dotted curves result from ilr = fi- (Compare the curves with 
the inset in Fig. 7 in (ill.) 
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FIG. 9: Same as Fig.O except for the black hole spin param- 
eter fir- 



star gets swallowed by the black hole more or less intact, 
albeit distorted by tidal forces. Clearly, this question has 
important consequences from the perspective of gravita- 
tional wave observations, but perhaps even more impor- 
tant are the ramifications for SGRBs. To launch such a 
burst requires the formation of an accretion disk around 
the black hole, which can occur only if the neutron star is 
disrupted prior to reaching an ISCO. To explore this issue 
quantitatively requires a dynamical simulation (compare 
d, [1^), and part of the motivation for the work pre- 
sented in this paper is the construction of suitable initial 
data for such calculations. In the meantime, however, 
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we may also use our quasiequilibrium models to obtain 
preliminary estimates. Specifically, we will use our nu- 
merical results to construct qualitative expressions that 
may be used to predict whether a BHNS binary of arbi- 
trary mass ratio and neutron star compaction encounters 
an ISCO before being tidally disrupted or not. 
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(see [4^)- This indicator takes the value unity for a 
spherical star and reaches zero at the formation of a cusp. 
We tabulate x as a function of the orbital angular veloc- 
ity and extrapolate to x = by using fitting polynomial 
functions to find the onset of tidal disruption. In Fig. [TOl 
we show an example of such extrapolations for sequences 
of neutron-star baryon rest mass Mb = 0.15 with mass 
ratios 1, 2, 3, and 5. Note that the horizontal axis in 
Fig. [TO]is the orbital angular velocity in polytropic units, 
Cl = i7i?poiy. Then we prepare the data of extrapolation 
for all models we compute. 

To obtain the fitting formula for the orbital angular 
velocities at the mass-shedding limit acceptable for all 
models we compute, we start with the qualitative New- 
tonian expression (pS]) . together with Kepler's third law 



(28) 



Combining Eqs. (|26p with we can eliminate dtid and 
find 



f^tid = 0.270 



^3/2 



(l + .) , (29) 



or equivalently 



FIG. 10: Extrapolation of sequences for the neutron-star 
baryon rest mass Mb = 0.15 to the mass-shedding point 
(x = 0)- Thick lines are sequences constructed using nu- 
merical data, and thin lines are extrapolated sequences. 

We start from investigating the binary separation (and 
the orbital angular velocity) at which tidal disruption 
of the neutron star by the black hole occurs. In New- 
tonian gravity and semi-relativistic approaches, simple 
equations may be introduced to fit the effective radius 
of a Roche lobe [2l|, [sl, [59| . Recently, Shibata and 
Uryii introduced a fitting equation for binaries composed 
of a nonspinning black hole and a corotating neutron star 
in general relativity 9]. In this section, we investigate 
our data for a nonspinning black hole and an irrotational 
neutron star. 

In order to obtain a fitting formula, we need to deter- 
mine the orbital angular velocity at the mass-shedding 
limit, i.e., the point that defines the onset of tidal dis- 
ruption. For this purpose we need to extrapolate our 
data, because it is impossible to compute sequences up 
to the tidal disruption point using a spectral methods 
code because of Gibbs phenomena. To extrapolate our 
results, we introduce a sensitive mass-shedding indica- 
tor, X, defined as the ratio between radial derivatives of 
the specific enthalpy h on the neutron star surface in the 
direction of the companion and in the polar direction. 



{d {In h)/dr) pole 



(27) 



f^tidMo = 0.270C3/2(1 + q)(l + (30) 



Here we have identified Mbh with M^^, Mns with 
-^^ADMO' Rns with the circumferential radius of a 
spherical neutron star Rq. Lastly, we determined the 
coefficient of 0.270 from our numerical results. Fitting 
our numerical values for the angular velocity at tidal dis- 
ruption to the form (j29|) resulted in a narrow range of 
coefficients between 0.266 and 0.273, with a mean value 
of 0.270. In Fig. [11] we show the results of the fitting of 
the mass-shedding limit by Eq. ((29|) . 

It may be of interest to compare Eq. (j29p with a sim- 
ilar expression of Shibata and Uryii [9| , who express the 
critical mass ratio q as 

(in the mass ratio g* is defined as the inverse of our 
definition). In terms of the quantities used in Eq. (j29|) . 
this becomes 



a 



tid,SU 



0.27 



^3/2 



jTfNS 



(32) 



which agrees with Eq. ((29)l for large q. The agreement of 
our result (irrotation for the neutron star) with that by 
Shibata and Uryii (corotation) confirms our prediction in 
the limit of extreme mass ratios 1341 . 
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FIG. 11: Fits of the mass-shedding limit by the analytic ex- 
pression (|29|) . The mass-shedding limit for each neutron-star 
mass and mass ratio is computed by the extrapolation of the 
numerical data. 



B. Innermost stable circular orbit 

We are also interested in the binary separation (and the 
corresponding orbital angular velocity) at which the min- 
imum of the binding energy appears, corresponding to 
the ISCO. In our numerical data we locate this point by 
fitting a second-order polynomial to three nearby points. 

We may construct a simple empirical fit that predicts 
the ISCO angular velocity rJisco for an arbitrary com- 
panion orbiting a black hole as follows. We search for ex- 
pressions with three free parameters that express ilisco 
as a function of the mass ratio q and the compaction 
C of the companion. We then fix the three parameters 
by matching to three known values of fiisco, namely 
(1) that of a test particle orbiting a Schwarzschild black 
hole, fiiscoMo = 6-3/2 ^fgj. ^ ^ ^-j^ that of an 
equal-mass binary black hole system as computed in [43| , 
^iscoMo = 0.1227 (for g = 1 and C = 0.5),_and finally 
(3) that of our BHNS configuration with M^^ = 0.15 
and C = 0.1452 and mass ratio q — 5, resulting in 
f^isco-^^o = 0.0854. A further requirement arises from 
the fact that for a test particle (with q = oo), the ex- 
pression should be independent of the companion's com- 
paction. We find a good fit to our remaining numerical 
data with the expression 

l^iscoAfo = 0.0680 (^1 - (l - 3.54Ci/3)^ , (33) 

as demonstrated in Fig. [121 We determined the expo- 
nent of q and C in Eq. (j33p empirically by requiring that 
the fitted curves lie near the data points for all models. 
Clearly the agreement is not perfect, but adequate for 
our purposes. 



FIG. 12: Fits of the minimum point of the binding energy 
curve by expression (|33p . 



C. Critical mass ratio 

Combining Eqs. ((29l) (or JSQl)) and ((33)l . we can now 
identify the critical binary parameters that separate 
those binaries that encounter an ISCO before reaching 
mass-shedding, and vice-versa. In Fig. [12] we show an 
example for M^^ = 0.15. The solid curve denotes the 
orbital angular velocity of the mass-shedding limit, and 
the long-dashed for the ISCO. As seen from Eqs. pS)) 
and (|33p . both of these curves depend on the mass ratio, 
but in different ways, which leads to the intersection of 
the two curves. An inspiraling binary evolves along hori- 
zontal lines towards increasing fj, starting at the left and 
moving toward the right, until reaching either the ISCO 
or the mass-shedding limit. After the binary reaches the 
ISCO for sufficiently large mass ratio, we cannot deter- 
mine the exact position at which the binary disrupts via 
quasiequilibrium calculations because it is in the dynam- 
ical plunge region, but we nevertheless include the mass- 
shedding limit for unstable quasiequilibrium sequences 
as the dotted curve in Fig. [T31 As shown in Fig. [131 the 
model with mass ratio q = 6 (dotted-dashed line) en- 
counters the ISCO, while that of g = 3 (dot-dot-dashed 
line) ends up at the mass-shedding limit. The intersec- 
tion between the mass-shedding and ISCO curve marks 
a critical point that separates the two distinct outcomes 
of the binary inspiral. 

Since expressions (I^H)) and also depend on the neu- 
tron star's compaction, the specific values of the critical 
point likewise depend on this compaction. In Fig. I14i we 
graph the mass-shedding and ISCO curves for a number 
of different neutron star masses which have a one-to-one 
correspondence to the compaction. The intersections, ap- 
pearing as a knee in this figure, mark the critical point 
for each compaction corresponding to a binary that en- 
counters the ISCO at the point of tidal disruption. This 
critical point also gives the maximum value of CIMq in 
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FIG. 13: An example of the boundary between the mass- 
shedding limit and the ISCO. We select the model of neutron 
star mass Afg^ = 0.15. The solid curve denotes the mass- 
shedding limit, and the long-dashed one the ISCO for each 
mass ratio as a function of the orbital angular velocity in 
poly tropic units. The dotted curve denotes the mass-shedding 
limit for unstable quasiequilibrium sequences. 

quasiequilibrium for each model. We note again that the 
equation of state is fixed here to be polytropic with adi- 
abatic index F = 2. This implies that the results may 
change, i.e., the coefficients in Eqs. and ([55]) may 
change, when we change the adiabatic index or the equa- 
tion of state itself. However, we find from Fig.[T3]that the 
orbital angular velocity at the ISCO in the form of fJA/o 
has a narrow range, 0.08 ^ fiMo % 0.09, for all mod- 
els we compute. This implies that the fitting formula 
P3p may hold approximately even for other equations of 
state. 

We tabulate these critical values for the neutron star 
compactions considered in this paper in TableU Whether 
or not these critical points mark the true separation be- 
tween the two outcomes will have to be verified by dy- 
namical simulations. We nevertheless expect that these 
values may provide some useful guidance. 

TABLE I: Critical value at which tidal disruption occurs at 
the ISCO for different neutron star models. 
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When we eliminate fiMo from Eqs. ([50)1 and ([55]) . we 
can draw a curve of the critical mass ratio which sep- 



FIG. 14: Endpoint of sequences for different neutron star 
masses. Left panel is the mass ratio as a function of the orbital 
angular velocity in polytropic units, while the right panel is 
normalized by the ADM mass at infinite orbital separation. 



arates BHNS binaries that encounter an ISCO before 
reaching mass-shedding, and vice-versa, as a function of 
the compaction of the neutron star. We show such a 
critical curve that separates those two regions in Fig. [151 
The solid line denotes the critical mass ratio for each 
compaction. If the mass ratio of a BHNS binary is larger 
than the critical one, the quasiequilibrium sequence ter- 
minates by encountering the ISCO, while if smaller, it 
ends at the mass-shedding limit of the neutron star. 
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FIG. 15: Critical mass ratio which separates BHNS binaries 
that encounter an ISCO before reaching mass-shedding and 
undergoing tidal disruption, as a function of the compaction 
of neutron star. 
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V. SUMMARY 

We have constructed new quasiequilibrium configura- 
tions of black hole-neutron star binaries in general rela- 
tivity. We have solved the Einstein constraint equations 
in the conformal thin-sandwich formalism coupled with 
the equations of relativistic hydrostationary equilibrium. 
In Paper I, we set the spin angular velocity parameter of 
the black hole equal to that of the orbital angular veloc- 
ity in order to produce a nonspinning black hole in the 
"leading-order" approximation ^41,], while in this paper 
we compute this parameter by requiring the quasilocal 
spin angular momentum of the black hole to be zero [i^ . 
We have also improved the formulation of the gravita- 
tional field equations and obtained more accurate results 
than in Paper I. 

As an indication of the improvements in these calcula- 
tions, a post-Newtonian analysis predicts smaller binary 
eccentricities for these new BHNS models than for those 
computed in Paper I ([s^], compare [6l|). In [6l|, Berti 
et at fit numerical results for the binding energy and 
angular momentum of binaries in circular orbit to post- 
Newtonian expressions for binaries that are not necessar- 
ily in circular orbit. Deviations between the the two ap- 
proaches then lead to non-zero eccentricities in the post- 
Newtonian expressions. These eccentricities are smaller 
for our new results than for those of Paper I. We also re- 
mark on another finding of (sij . namely that for a given 
neutron star mass M^dm o ^^"^ ^ given value of flMo, the 
eccentricities in BHNS models, though small, are found 
to be larger than in binary neutron star models [63 |. 
This suggests a larger deviation from quasiequilibrium 
for BHNS binaries than binary neutron stars. But, for 
BHNS binaries with a mass ratio of q = 5, these param- 
eters correspond to a larger binary separation than for 
binary neutron stars with a mass ratio oi q = 1 (compare 
Eq. ((28)) ). For similar numerical resources, this larger bi- 
nary separation leads to a larger numerical error, which 
may explain the larger eccentricity found by |6l| , at least 
in part. 

In addition to recomputing several sequences we pre- 



sented in Paper I, we have constructed sequences for a 
wider range of neutron star masses and binary mass ra- 
tios, employing a F = 2 polytropic neutron-star equation 
of state throughout. We computed several constant-mass 
sequences, for various mass ratios and neutron star com- 
pactions, and searched for the appearance of a cusp at 
the neutron star surface - indicating the onset of tidal 
disruption - and turning points on the binding energy 
and angular momentum curves - identifying the ISCO. 
We also included some qualitative fits that allow for a 
simple prediction of those binary parameters separating 
these two different outcomes of binary coalescence. Un- 
like in our earlier findings, we found simultaneous turning 
points along the binding energy and angular momentum 
quasiequilibrium curves. 
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APPENDIX A: TABLES OF SEQUENCES 

We summarize our results in Tables |lT] - EH In these 
tables, we tabulate the coordinate orbital separation be- 
tween the center of the excised surface of the black hole 
and the maximum of the neutron-star baryon rest mass 
d, orbital angular velocity Q, spin angular velocity of the 
black hole fir, binary binding energy iJb, total angular 
momentum J, decrease in the maximum density param- 
eter (5(7niax = (gmax " gmax,o)/9max,o from the Spherical 
value (/max.Oj minimum of the mass-shedding indicator 
Xinin, and fractional difference 5M between the ADM 
mass Madm and the Komar mass MKomar along a se- 
quence. Here (?max = (-F'/po)max- Recall that (virial) 
equilibrium requires Af Komar — Mabm- 
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TABLE III: Same as Table |TT] but for the neutron-star baryon rest mass A'/f^ = 0.16. The ADM mass and the isotropic 
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TABLE IV: Same as Table |II] but for the neutron-star baryon rest mass M^^ = 0.15. The ADM mass and the isotropic 
coordinate radius of the neutron star in isolation are Madm.o = 0.1395 and fo — 0.8152 (k = 1). The neutron star compaction 
is C = 0.1452. 



d/Mo QMo 



Mass 

QrMo 



ratio: 
Eb/Mo 



f2 



= 1 



J/M^ Sqn 



Xmin 



5M 



15.28 
13.10 
10.91 
9.82 
9.17 
8.73 
8.51 
8.40 



0.0154 
0.0191 
0.0246 
0.0285 
0.0313 
0.0335 
0.0347 
0.0353 



0.0145 
0.0179 
0.0227 
0.0261 
0.0285 
0.0303 
0.0313 
0.0318 



-7.30(-3) 
-8.34(-3) 
-9.70(-3) 
-1.06(-2) 
-l.ll(-2) 
-1.15(-2) 
-1.17(-2) 
-1.19(-2) 



1.109 
1.049 
0.986 
0.954 
0.934 
0.922 
0.916 
0.913 



-5.93 
-9.17 
-1.86 
-3.20 
-4.83 
-6.66 
-8.04 
-9.05 



0.961 1.92(-5) 
0.934 1.20(-5) 
0.871 1.27(-6) 
0.803 2.90(-6) 
0.731 1.19(-5) 
0.648 2.77(-5) 
0.574 4.80(-5) 
0.514 7.56(-5) 



/■NS 
^ADM,( 



d/Mo nMo 



Mass ratio: M§^/Mj 



Xmin 



5M 



20.37 
11.64 
8.73 
8.00 
7.28 
6.84 
6.55 
6.48 



0.0102 
0.0225 
0.0334 
0.0376 
0.0427 
0.0464 
0.0491 
0.0498 



0.00989 
0.0213 
0.0311 
0.0348 
0.0392 
0.0423 
0.0446 
0.0452 



-4.99(-3) 
-8.21(-3) 
-1.03(-2) 
-1.10(-2) 
-1.17(-2) 
-1.22(-2) 
-1.26(-2) 
-1.27(-2) 



1.104 
0.894 
0.818 
0.799 
0.781 
0.770 
0.763 
0.761 



-7.31 
-2.48 
-5.02 
-6.56 
-9.39 
-1.25 
-1.58 
-1.69 



0.993 
0.952 
0.870 
0.822 
0.742 
0.661 
0.565 
0.525 



3.67(-5) 
5.67(-5) 
6.10(-5) 
5.38(-5) 
5.84(-5) 
8.16(-5) 
1.24(-4) 
1.48 (-4) 



Mass ratio: M^^/M^^ma 



d/Mo 


Q.Mo 


fir Mo 


Eb/Mo 


J/Ml 


dqm 


IX Xmin 5M 


19.65 


0.0107 


0.0105 


-4.29(-3) 


0.918 


-1.22 


-3) 0.996 1.36(-4) 


13.10 


0.0191 


0.0184 


-6.22(-3) 


0.786 


-2.70 


-3) 0.982 4.32(-5) 


9.82 


0.0284 


0.0270 


-7.94(-3) 


0.715 


-5.11 


-3) 0.953 4.76(-5) 


7.64 


0.0399 


0.0375 


-9.61(-3) 


0.666 


-9.14 


-3) 0.888 1.01(-4) 


6.55 


0.0490 


0.0455 


-1.06(-2) 


0.643 


-1.35 


-2) 0.808 1.09(-4) 


6.00 


0.0549 


0.0506 


-1.12(-2) 


0.633 


-1.76 


-2) 0.734 1.19(-4) 


5.46 


0.0620 


0.0567 


-1.18(-2) 


0.624 


-2.53 


-2) 0.595 1.80(-4) 


5.35 


0.0636 


0.0581 


-1.19(-2) 


0.622 


-2.78 


-2) 0.537 2.17(-4) 






Mass ratio: M 


BH/,rNS _ 
rr /^"ADM,0 — 


5 


d/Mo 


Q,Mo 


QrMo 


E^/Mo 


J/Ml 


dqm 


IX Xmin 5M 



18.19 0.0120 0.0118 -3.32(-3) 0.6605 -2.19 

11.64 0.0224 0.0218 -4.98(-3) 0.5601 -4.58 

7.28 0.0425 0.0406 -7.22(-3) 0.4882 -1.29 

5.46 0.0618 0.0580 -8.60(-3) 0.4624 -2.53 

4.59 0.0769 0.0713 -9.12(-3) 0.4554 -3.92 

4.37 0.0817 0.0755 -9.177(-3) 0.45484 -4.51 

t 4.22 0.0851 0.0784 -9.192(-3) 0.45485 -5.02 

4.12 0.0879 0.0808 -9.185(-3) 0.45500 -5.50 



0.999 
0.989 
0.945 
0.853 
0.730 
0.674 
0.625 
0.570 



2.48(-4) 
1.89(-4) 
8.70(-5) 
1.55(-5) 
1.89(-5) 
2.30(-5) 
3.25(-5) 
4.01(-5) 
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TABLE V: Same as Table |II] but for the neutron-star baryon rest mass Mf^ = 0.14. The ADM mass and the isotropic 
coordinate radius of the neutron star in isolation are M^dm.o = 0.1310 and fo — 0.8556 («: = 1). The neutron star compaction 
is C = 0.1321. 







Mass ratio: 


Mj 


BH/,^NS _ 
rr /-'"ADM.O — 


3 








d/Mo 


mio 









J/Mi 




Xmin 


SM 


13.95 


0.0174 


0.0167 


-5.91(- 


■3) 


0.805 - 


■2.17(- 


■3) 


0.978 


6.26(- 


■6) 


10.46 


0.0260 


0.0249 


-7.56(- 


■3) 


0.729 - 


■4.18(- 


■3) 


0.943 


7.02(- 


■5) 


8.14 


0.0367 


0.0346 


-9.18(- 


■3) 


0.678 - 


■7.78(- 


■3) 


0.865 


1.09(- 


-4) 


6.98 


0.0451 


0.0420 


-1.02(- 


■2) 


0.653 - 


■1.23(- 


■2) 


0.766 


1.18(- 


-4) 


6.63 


0.0482 


0.0448 


-1.05(- 


■2) 


0.646 - 


■1.48(- 


■2) 


0.714 


1.28(- 


-4) 


6.39 


0.0505 


0.0468 


-1.07(- 


■2) 


0.642 - 


■1.70(- 


■2) 


0.668 


1.43(- 


-4) 


6.16 


0.0530 


0.0490 


-1.10(- 


■2) 


0.638 - 


■2.01(- 


■2) 


0.604 


1.70(- 


-4) 


5.99 


0.0551 


0.0508 


-1.12(- 


-2) 


0.635 - 


■2.32(- 


■2) 


0.530 


2.06(- 


-4) 






Mass ratio: 


Mi 


rr /-'"ADM.O — 


5 









d/Mo QA^o nrMo Eb/Mo J/Mp Sqm^^ x^in 5M 



10.85 


0.0247 


0.0240 


-5.36(- 


■3) 


0.547 


-4.95(- 


■3) 


0.979 


6.41(- 


■5) 


7.75 


0.0391 


0.0374 


-6.96(- 


■3) 


0.496 


-1.08(- 


■2) 


0.934 


3.33(- 


-6) 


6.20 


0.0524 


0.0496 


-8.08(- 


■3) 


0.473 


-1.85(- 


■2) 


0.859 


8.52(- 


■5) 


5.43 


0.0622 


0.0584 


-8.64(- 


■3) 


0.463 


-2.61(- 


■2) 


0.777 


1.12(- 


-4) 


5.04 


0.0683 


0.0639 


-8.89(- 


■3) 


0.460 


-3.23(- 


■2) 


0.707 


1.25(- 


-4) 


4.73 


0.0740 


0.0688 


-9.06(- 


■3) 


0.458 


-4.00(- 


■2) 


0.619 


1.47(- 


-4) 


4.58 


0.0771 


0.0716 


-9.12(- 


■3) 


0.4572 


-4.54(- 


■2) 


0.546 


1.65(- 


-4) 


4.54 


0.0780 


0.0723 


-9.14(- 


■3) 


0.4570 


-4.72(- 


■2) 


0.516 


1.77(- 


-4) 



TABLE VI: Same as Table |II] but for the neutron-star baryon rest mass Af^^ = 0.13. The ADM mass and the isotropic 
coordinate radius of the neutron star in isolation are Madm,o = 0.1223 and fo = 0.8923 (k = 1). The neutron star compaction 
is C = 0.1201. 



Mass ratio: M^^/M^^^.o = 3 

d/Mo ^Mq nrMo gb/Mp J/M^ ggmax Xmin SM 



14.93 


0.0158 


0.0154 


-5.58(- 


■3) 


0.825 


-1.79(- 


■3) 


0.974 


1.39(- 


■5) 


11.20 


0.0237 


0.0227 


-7.15(- 


■3) 


0.747 


-3.51(- 


■3) 


0.933 


7.52(- 


■5) 


8.71 


0.0334 


0.0316 


-8.71(- 


■3) 


0.691 


-6.94(- 


■3) 


0.841 


1.03(- 


-4) 


8.09 


0.0370 


0.0348 


-9.19(- 


■3) 


0.678 


-8.82(- 


■3) 


0.792 


1.05(- 


-4) 


7.47 


0.0412 


0.0386 


-9.71(- 


■3) 


0.665 


-1.19(- 


■2) 


0.717 


1.16(- 


-4) 


6.97 


0.0451 


0.0421 


-1.02(- 


■2) 


0.655 


-1.63(- 


■2) 


0.617 


1.45(- 


-4) 


6.85 


0.0462 


0.0431 


-1.03(- 


■2) 


0.652 


-1.79(- 


■2) 


0.580 


1.60(- 


-4) 


6.72 


0.0473 


0.0441 


-1.04(- 


■2) 


0.650 


-1.97(- 


■2) 


0.531 


1.79(- 


-4) 



Mass ratio: M™/Madm,o = 5 



d/Mo i^Mo rirMp E^/Mo J/Mj Sq^^y, Xmin SM 



11.62 


0.0225 


0.0219 


-5.06(- 


■3) 


0.560 


-4.13(- 


■3) 


0.975 


6.73(- 


■5) 


8.30 


0.0356 


0.0342 


-6.66(- 


■3) 


0.505 


-9.10(- 


■3) 


0.923 


5.36(- 


■5) 


7.47 


0.0410 


0.0392 


-7.19(- 


■3) 


0.492 


-1.18(- 


■2) 


0.890 


9.24(- 


■5) 


6.64 


0.0479 


0.0456 


-7.76(- 


■3) 


0.480 


-1.61(- 


■2) 


0.835 


1.26(- 


-4) 


5.81 


0.0570 


0.0538 


-8.36(- 


■3) 


0.469 


-2.36(- 


■2) 


0.735 


1.60(- 


-4) 


5.40 


0.0627 


0.0589 


-8.65(- 


■3) 


0.465 


-3.06(- 


■2) 


0.645 


1.90(- 


-4) 


5.23 


0.0652 


0.0611 


-8.76(- 


■3) 


0.463 


-3.45(- 


■2) 


0.590 


2.09(- 


-4) 


5.11 


0.0672 


0.0629 


-8.84(- 


■3) 


0.462 


-3.81(- 


■2) 


0.532 


2.24(- 


-4) 



